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We have applied the transformation of the slow light equations to Liouville theory that we devel- 
oped in our previous work, to study the influence of relaxation on the soliton dynamics. We solved 
the problem of the soliton dynamics in the presence of relaxation and found that the spontaneous 
emission from the upper atomic level is strongly suppressed. Our solution proves that the spatial 
shape of the soliton is well preserved even if the relaxation time is much shorter than the soliton 
time length. This fact is of great importance for applications of the slow-light soliton concept in 
optical information processing. We also demonstrate that the relaxation plays a role of resistance 
to the soliton motion and slows the soliton down even if the controlling field is constant. 
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The development of modern methods of optical sig- 
nal manipulation and control opens wide perspectives for 
application of the light in classical and quantum compu- 
tation. One of such methods of optical signal process- 
ing is based on the effect of electromagnetically induced 
transparency (EIT) [l| , which allows for slowing the light 
down by many orders of magnitude and even bringing it 
to a complete halt 0,0, 0]. The major advantage of this 
methods is that the velocity of the optical signal is ef- 
fectively controlled by an auxiliary laser [1, @, 0] . In the 
linear regime of operation when the intensity of the con- 
trolling field is significantly larger than the intensity of 
the signal, there are some important constraints imposed 
on the parameters of the signal [l| . These constraints re- 
sult from the presence of strong optical relaxation in the 
medium. The medium is typically a gas of alkali atoms 
whose electronic structure relevant to the EIT effect can 
be schematically described by the three-level A model 
(see Fig. [T]). The relaxation results from spontaneous 
transitions of the atoms from the excited upper energy 
level 1 3). Usually, such transitions may occur not only to 
the levels |1),|2) but also to other lower levels, which are 
not included into consideration for the sake of simplicity. 
In any case, the relaxation destroys optical coherence of 



the signal and therefore must be accounted for. 

The atom-field interaction in the system with EIT is 
substantially nonlinear. In our previous works we pro- 
vided a nonlinear solution of the EIT problem called the 
slow- light soliton We also demonstrated that the 
dynamics of slow-light solitons strongly depends on the 
form of the controlling field and the solitons can be effec- 
tively manipulated concerning their space-time dynam- 
ics [i| [HI. In our present work we study the dynam- 
ics of the solitons in the presence of the relaxation. We 
note that in the framework of the nonlinear theory there 
might exist different approaches to effectively suppress 
the influence of relaxation. Here, we demonstrate that 
appropriately choosing the parameters of the controlling 
field we can reduce the influence of the relaxation on the 
slow-light solitons by several orders of magnitude. Our 
results are found to be in a good agreement with the 
experiments on EIT. 

Within the slowly varying amplitude and phase ap- 
proximation (SVEPA) the Maxwell equations for the 
Rabi- frequencies are reduced to the the well-known equa- 
tions governing the dynamics of the atom-field system, 
viz. 
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Here £ = zj 'c,t — t — z / c. The Schrodinger equation for 
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|3)= 3P 3/2 ,F = 2,M F = 0) 




|l) = |3S V2 ,F = l,M F = -l) 
FIG. 1: The A-scheme for sodium atoms. 

the amplitudes V'1.2,3 of atomic wave function reads 



2 n l ^3; 



(2) 



(i A + 1) V-3 + ^(O^l + Qfc^)- 



By excluding the amplitudes of the lower levels the sys- 
tem of Eqs.©, © can be transformed to the following 
form 



5 T |^3| 2 =-7l^3| 2 -^ c (|fia| 2 

d T ifi 3 = -A + 2t/o ^ |a (|r> a | 2 d^ a + 



(3) 
(4) 

H^l 2 )- (5) 
n b \ 2 d cVb ). (6) 



Here, ip a ,b,3 are respectively the phases of the fields fl a b 
and the amplitude of the excited state ^3. Notice that 
the equations Eqs. ©,© are the wave equations for the 
fields in curvilinear coordinates with the metric depen- 
dent on the amplitude ^3. 

In the typical dynamics scenario, the atoms are ini- 
tially prepared on the lowest energy level |1), the probe 
field is absent while the controlling field is constant. No- 
tice that the state |1) is a dark state for the controlling 
field. Thus, the initial state of the system reads 



n„(o,c) = 0, n 6 (o,o = n , |V.t(o,0) 
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(7) 



In fact, this is the simplest solution of the equations 
Eqs.©,©,©,©. 



I. 



RELAXATION 



From the results of linear and nonlinear theories of 
electro mag netically induced transparency (EIT) @, [ll|, 
[H, [HI, ll4L [l5| we infered a physically plausible assump- 
tion that the population of the excited state |3) is propor- 
tional to the amplitude of the field in the probe channel, 



which we denote as fl a ■ In our previous works we demon- 
strated that a slow-light soliton in the system without 



relaxation is obtained by taking ^3 = 



jQ a . Here, 



2|A-A| 

A is an arbitrary parameter limited from below by the 
condition | 1 < 1- ln order to take into account the re- 
laxation we propose a more general relation between the 
amplitude of atomic upper level |3) and the amplitude 
il a of the slow-light soliton. The population of the up- 
per level although very small should remain nonzero in 
order to preserve necessary coupling between the probe 
and controlling field. At the same time, the relaxation 
should manifest itself in the form of some effective damp- 
ing of the probe field. Based on such consideration we 
found the following relation between the amplitudes: 



o a ( T ) 



fp3 



1 



2IA- A| 



rO, 



(8) 



where a(r) is a first order correction to the exact slow- 
light soliton solution. It will be shown below that a(r) is 
always negative. Let us introduce convenient notations 



= e -p 



= 77, p = p + a, k = 



IA 



The first three equations of the system Eq.(??) take then 
the following form: 



d T a d^p + d(; T p = —k e 2p 
d( T r] + d T p d^T] = k e^ 2 ? 77, 
Ak(d T + 7) e~ 2 P = -d c (rf + e~ 2 P) 



(9) 
(10) 
(11) 



The last equation for the phase tp 3 can easily be inte- 
grated after we solved the first three. 

Assuming that the effective relaxation described by a 
varies slowly in time r and therefore neglecting the first 
term in Eq.© we can find a solution of the system of 
equations Eqs.©, (fTUf . (fTTj) . We will provide a rigorous 
description of this approximation after we construct the 
solution in the next section. 



II. SOLUTION 

We assume that d T a is small as compared to the r.h.s. 
of the equation Eq. © , and neglect the first term in the 
l.h.s of this equation. The equation then transforms into 
the well known Liouville equation, viz. 



d C rp 



-k e 



-2/5 



whose general solution is readily available: 



log 



le 2a d c A + (C)d T A_( T ) 



(12) 



(13) 



Here A + ((), A_ (t) are arbitrary functions. To obtain 
the solution of the whole system including Eqs. (fTT)|) and 
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(fTTjl we specify these arbitrary functions as follows: 

A+(C) =-exp[-8£ fcCL 
A_(r) = exp 2e / Jt^ dr 
i] = —2pd T p + 2d T p — 2pd T a 



d T a(r) = — 



7/2 

p(r)2 + l • 



(14) 
(15) 
(16) 

(17) 



p(r) is an arbitrary function describing the controlling 
field n(r), while A = ie - We note from Eq. frT)) that 
the correction a to exact solution without relaxation ob- 
tained in our previous papers vanishes for 7 = as ex- 
pected. 

So, the fields read 



O»(r,0 



sech(^), (18) 

28 T p(r)-7p(r) 



with the phase 



</3 = -4fce (C — Co) +£o : 



2c(r) 



(19) 



We chose a(0) = 0. From the solution for the phase of 
slow-light soliton Eq.(|19|) we obtain its velocity: 



(20) 



4fcp(r) 



From equation Eq. (|22|) we can find the spatial window of 
validity of our approximation with respect to the center 
of slow- light soliton A£(t) = C — Cc(t) where our soliton 
provides correct description of the pulse shape. It reads 



1 (\f 

Isinhf-Sfceo AC) I < — e 2a . 

7 



(25) 



At the initial moment of time, for 7 < 16eo we find 
ke A( « In (2a). We will denote the full-width at 

half-maximum of slow-light soliton by w s w 0.66/(fc|en|). 
To keep our solution valid at least within the full-width 
at half-maximum of the soliton, i.e. w s = 2ACo, the 
parameter eo should obey the following condition 



£0 > 0.7 7 



(26) 



Notice that a is a negative monotonically decreasing 
function of r. Therefore, the validity window closes down 
with the retarded time r. However, the group velocity of 
the slow-light soliton approaches zero as r increases and 
the soliton slows down to full stop (see Eq. (|20|) ). The 
distance that the soliton travels until full stop is 



£ = C c (oo)-Co<-^. 



(27) 



This formula clearly shows that the maximum distance 
that our soliton can travel in the medium is limited by 
the magnitude of the relaxation constant. As stronger 
the relaxation is, as smaller the distance that the soliton 
can propagate in the medium. 



III. VALIDITY OF THE APPROXIMATION 

Now we will analyze the validity of our approximation 
of small d T a. The first term of Eq.© can be neglected 
if it is much smaller than ke~ 2 ?. Hence we find the 
following condition on a: 



\d T a\ < 



8soe 



2a 



(p 2 + 1) sinh(2y>) 



(21) 



This condition can be simplified with the help of Eq. (fT7|) 
and we obtain 



7 



16e 



< 



sinh(2<^)| 



(22) 



It is always fulfilled at the maximum of slow-light soliton, 
because ip = there. The position of the maximum of 
the soliton is given by the following function of retarded 
time r: 



1 _ e 2a(r) 



(23) 



where we chose Cc(0) = Co- Hence, the phase of the 
soliton can be rewritten as 



^=-4fce (C-Cc(r)). 



(24) 



IV. COMPARISON WITH EXPERIMENTS 

We discussed in our previous papers that the control- 
ling field generated by an auxiliary laser on entrance into 
the medium supports the propagation of the slow-light 
soliton. An appropriately chosen modulation of the con- 
trolling field ensures the nonlinear coupling between the 
soliton and atoms, while the nonlinear coupling in its 
turn preserves the secant shape of the soliton. However, 
in the ideal case when the relaxation is absent, the mod- 
ulation of the controlling field has the same time length 
as the signal in the channel a. After the soliton has 
been generated, the auxiliary laser continues emitting a 
constant laser beam, which supports the propagation of 
the soliton with a constant velocity. We demonstrated 
that this ideal situation when the amplitude of the soli- 
ton remains constant corresponds to a constant function 
p(r) = Po- If the relaxation is not zero, the amplitude 
and the velocity of the soliton will be decaying in time r 
according to Eq. (fT8|) . The asymptotic value of the back- 
ground field past the modulation supporting the soliton 
can be found from r2b(r, C) by taking the limit C — * 00 at 
the zero moment of time: 



r» = 



(2|g |-7)po 



(28) 
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FIG. 2: Slow-light soliton decay due to the relaxation. The 
solid line corresponds to a reference soliton propagating in the 
system without relaxation, and the dashed line represents the 
pulse in the presence of relaxation. 



Assuming that 7 3> ilo we find 



Po 




(29) 



In the absence of the relaxation, the time length of the 
soliton pulse on entrance into the medium is determined 
by the phase <po(t) = eo*/(Po + !)■ Hence we obtain a 
second condition on two arbitrary parameters po and Sq 
in terms of well-defined experimental parameters f^o and 

t n • 



sech 



£0 



Pi 



1 



t P /2 



0.5 



(30) 



Choosing |eo| = 7 we verify that the conditions Eqs. 
(j28|) . (|30p are well satisfied for experimentally feasible 



values of t p ~ lfis and f^o ~ 10 MHz. 

Now we will compare our solution with the experimen- 
tal results reported for sodium atoms 0]. The param- 
eters have the following values: 7 = 6.3 x 10 7 rads -1 , 
t p = 2.5fis, and Q = O.567. We solve Eqs.((25D, (HD 
and find p ~ 18.4, e — 5.77. We can also calculate the 
reduction in the strength of optical relaxation influenc- 
ing the dynamics of the slow-light soliton. According to 
our solution Eq. |[T5|) , the amplitude of the soliton decays 
with the rate 2a. Therefore, we can introduce an effective 
relaxation constant defined as follows (see Eq.(fl7|0 



7 



7 



Pi 



1 



(31) 



In the case considered above, 7* ~ 7/340. We note that 
the effective value of the relaxation constant 7 is signif- 



icantly lower for the soliton than for an arbitrary pulse. 
So, the spontaneous emission is greatly suppressed for 
soliton pulses due to the nonlinear interaction with the 
medium. The effective optical relaxation time r* el is 
larger than the pulse length and approximately equals 
2.2 t p . For the pulse delay At — 7.05 [is reported by 
Hau we calculated the decay of the soliton pulse am- 
plitude at the maximum and compared it with a reference 
pulse Q* propagating in the system without relaxation. 
We obtained 



n a {AT + t p /2Xc) 



0.2 



The reference pulse f2* is modeled traveling under the 
same conditions in the system without relaxation. Wc 
added a half of the pulse time width, because it started 
interacting with the medium approximately t p /2 mi- 
croseconds earlier than the maximum entered inside 
Fig. [Q This result agrees very well with the measure- 
ments reported by Hau. Wc also note that the distance 
that the slow-light soliton propagated in the medium dur- 
ing the time Ar-\-t p is approximately equal 200/im, which 
is of order of the atomic cloud in the experiment [2|. 
We need to emphasize that in the presence of the relax- 
ation the velocity of the soliton is not a constant any 
more like in the ideal case (see Eq. (|2"0)) . We calculated 
that the average value of the velocity is approximately 
22m I s compared to experimentally measured 32to/s for 
the Gaussian pulse. 

A more important feature of the dynamics is in the 
fact that the soliton spatial form does not depend on the 
relaxation at all. The soliton shape remains unchanged 
along the propagation inside the medium (see Fig. [IJ. 
This fact plays a crucial role for storage and retrieval of 
optical information in potential applications. 



V. DISCUSSION 

The results of our work remain valid beyond the con- 
straints of the transparency window defined in the linear 
theory. We demonstrated that due to strong nonlinear 
interaction between the probe and controlling fields it 
is possible to preserve the spacial shape of the optical 
signals even in the presence of strong optical relaxation. 
The comparison of our theory and experimental results 
0, [1[ shows a very good agreement. We provided rigor- 
ous analytical estimates for the largest distance that the 
slow-light soliton can propagate in the medium with re- 
laxation and also described the dependence of the soliton 
velocity on the relaxation constant 7. 
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